arXiv:1503.05745vl [math.AP] 19 Mar 2015 


A KINETIC REACTION MODEL: DECAY TO EQUILIBRIUM 
AND MACROSCOPIC LIMIT 

LUKAS NEUMANN, CHRISTIAN SCHMEISER 

Abstract. We propose a kinetic relaxation-model to describe a generation-recom¬ 
bination reaction of two species. The decay to equilibrium is studied by two recent 
methods for proving hypocoercivity of the linearized equations. Exponen¬ 

tial decay of small perturbations can be shown for the full nonlinear problem. 

The macroscopic/fast-reaction limit is derived rigorously employing entropy de¬ 
cay, resulting in a nonlinear diffusion equation for the difference of the position 
densities. 
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1. Introduction 

1.1. The model. We consider the system 

dtf + V - = Xi{v) - Pgf , 

dtg + v-Va:g = X2{v)-pfg, 

where / and g depend on position x G T^, the three dimensional torus of volume one, 
on velocity u G and on time f > 0. They represent the phase space densities of 
chemical reactants A and B, which are produced (with nonnegative velocity profiles 
Xi and X2, respectively) by the decomposition of a substance C. The density of the 
substance C is not subject of our study and is assumed to be fixed. On the other 
hand the substances A and B can recombine to form C and thus be eliminated from 
our system. Similar models have been used for generation and recombination of 
electron-hole pairs in semiconductors EE]. 

The probability of the reaction is depending on the position density 

Ph{x,t) := / h{x,v,t) du, 

7r3 

of the reaction partner. We consider the system ([1]) subject to initial conditions 

(2) f{x,v,0) = fo{x,v), g{x,v,0) = go{x,v), 
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with nonnegative initial data /o and g^. Moreover, since we want to describe the 
reaction lA+lBolC, we reqnire that ~ X 2 )dv = 0. We assnme that the 

system has been nondimensionalized and that Xi and X 2 satisfy 


(i + jvnxjeL\R^)nL^(R^), Xj>0, 

(3) / Xjdv = 1 , / vXjdv = 0 , 

Jr3 

3C, 6*>0: VaGM,a;G§^,(5>0: / Xjdv<C6^, j = l,2. 

J |a+f;-cj|<5 

The last line will be needed for an averaging lemma with the weight l/xj- The 
largest valne to be expected for the exponent is 0 = 1, which is achieved, e.g., for 
Gaussian distributions, the prototypical examples for the Xjj but more generally 
also for Xj{'^) ^ c(l + with k > 1. 

Note that, at least formally, the mass difference is conserved: 


d 

dt 



(/ - a) du dx 


0 , 


as can be seen by subtraction of the two equations and subsequent integration. 
This is to be expected since by the reaction molecules of A and B are created and 
destroyed pairwise. We introduce the unique poo > 0, such that 

(4) [ [ if 0 - go)dvdx = |T^| (p^ - —^ , 

Jt3 \ Poo / 


and expect convergence as f —?■ oo of solutions of ([I]), ([2]) to the steady state 
foe (,X,v) = PooXl (v) , goe{x,v) = —X2 (v) . 

poo 


This is supported by the decay properties of the entropy functional 

(5) H{f,g) = [ f (in— - l] + g (in^^ - l] dndx , 

J V PooXl J \ X2 J 


which decreases as long as (/,(?) is different from (p(x)xi(u), X 2 (u)/p(x)) for some 
p{x): 

(6) = j j /Ua2-/9')ln^dt,'d„di<0. 

Among these functions {foo^goo) is the only solution of ([T]) with the same mass 
difference as the initial data. 
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Spectral stability of the equilibrium will be investigated by linearization: 


(7) 


dtf + v ■ Vx/ = -PooXlPg -/ , 

poo 

dtg + v ■ = - X 2 pf - poog , 

Poo 


where for simplicity the perturbations have again been denoted by / and g, now 
satisfying 

(8) [ [ if - g)dv dx= [ (pf - pg)dx = 0 . 

JT3 JR3 Jtts 


Rigorous results for kinetic equations for chemically reacting species with nonlin¬ 
ear reaction models are scarce in the literature. An example is an existence result 
[13] for a model with quadratic reaction terms under rather weak natural assump¬ 
tions on the initial data, where stability is based on entropy decay. Since existence 
of solutions is not the main focus of this work, we shall make rather strong assump¬ 
tions on the data, with the consequence that the existence and uniqueness proof in 
the following section is based on weighted L°° estimates and rather straightforward. 

The analysis of the decay to equilibrium is complicated by the fact that the entropy 
dissipation (I6|) vanishes not only for the equilibrium, but on a larger set of local 
equilibria. If exponential decay to equilibrium can still be proven, the system is called 
hypocoercive [15]. The authors have been involved in the development of two different 
abstract procedures for proving hypocoercivity for linear equations [H 112], both 
based on the construction of suitable Lyapunov functionals (or modified entropies), 
whose dissipation functionals have appropriate coercivity properties. The method of 
[H] is based on a slightly tilted, weighted L^-norm, while [12] works in a setting 
and can be extended to higher regularity. In Section 2 we show that both methods 
are applicable to a linearized version of ([T]). Since the estimates of the existence 
result already provide neutral stability of the equilibrium, the decay results can be 
extended to a local asymptotic stability resnlt with exponential convergence for the 
full nonlinear model. The decay rates proven by both methods can, in principle, 
be compnted explicitly. Complete formulas would however be rather complicated, 
whence we did not attempt a comparison. An essential difference between the 
methods is the weaker assnmptions on initial data in [S]. On the other hand, the 
method of [12] has the potential to provide strong convergence properties inclnding 
derivatives. 

In Section [T75I the macroscopic/fast-reaction limit is carried ont formally, leading 
to a nonlinear diffusion equation for the difference of the position densities of the 
reactants. Similar results have been derived for reaction-diffusion systems [21E] 
and for coagulation-fragmentation models A rigorous justihcation of the limit 
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is the subject of Section 3. It is based on an analysis of the entropy dissipation 
functional ([H]) and adapts the procedure of [2] , where compactness is obtained from 
an averaging lemma in weighted L^-spaces. We prove a slightly generalized version 
compared to [2] . 

We note that with the torus we chose the simplest geometric setting. Natural 
modihcations include bounded domains with specular reflection boundary conditions 
or whole space problems with conhning potentials. We conjecture that our results 
can be extended to these situations, however with considerably more technical effort 
for the latter (see, e.g., the hypocoercivity results with conhning potentials in |H]). 


1.2. Existence of solutions. The entropy decay relation (E]) would suggest an 
existence result for initial data with bounded entropy. Such a result for a similar 
problem has been proven in [13]. The main ingredients are entropy inequalities, weak 
compactness and velocity averaging. These ideas might be transferable to our 
situation. However, for our purposes we need more information on the solutions. 
Under stronger assumptions on the initial data, a global existence result can be 
proved easily. 

Theorem 1. Let (Ej) hold and let there be positive constants 71 < poo 72 such 
that the initial data fo, go E L°°(T'^ x M^) satisfy 

(Poo - 7i)Xi < /o < {poo + l 2 )xi and < 9o < ■^;^X 2 ■ 

Then the initial value problem o. o has a unique global mild solution (/, g) G 
U([0, cxo), X satisfying 

{poo - 7i)Xi{v) < f{x,v,t) < {poo + 72 )Xi{v) , {x,v,t) e xM^ X [ 0 , 00), 

and 


^X2(u) < g{x,v,t) < ^^X2{v ), 


{x, V, f) G X X [0, cx)). 


Proof. The mild formulation of the equation for / is given by 


f{x,v,t) =fo{x-vt)exp 


Pg{x + v{s — f), s)ds 


+ Xi{'a) exp pg{x + v{t - t), T)dT^ ds , 


and an analogous equation holds for g. It is easily seen that the set of {f,g) de- 
hned by the estimates of the theorem is mapped into itself by the right hand sides. 
This provides the a priori estimate needed for the continuation of a local solution 
constructed by Picard iteration. □ 
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1.3. Formal macroscopic limit. In this section we formally derive a macroscopic 
limit of ([T]). The limit will be made rigorous in Section [31 Since by (|3]) the mean 
velocities of the equilibrium distributions vanish, we adopt a diffusive (or parabolic) 
scaling t ^ t/e^ and x ^ x/e and derive 

(?dtf+ = Xi{v)- pgf 

e^dtg + ev • W= X2(v) - pjg . 

We substitute the ansatz 


/ = r + e/' + 0{e^) and g = + eg^ + 0{e^) . 

Balancing the leading order terms gives Pgof^ = Xi and Pfog^ = X 2 - This is equiva¬ 
lent to the existence of p{x, t) such that 

f{x,v,t)=p{x,t)xi{v) and g^{x,v,t) =-^^-^X 2 {v). 

Now we balance the hrst order terms in e and derive 

P 

V ■ = -pfig° - pg^ . 

Due to (|3]) the solvability condition v ■ Va;(/° — g^)dv = 0 is satished, and we 
obtain 


f = -pXiv ■ VxP + p^Xi , 


g^ = -^X2V ■ Vxp - 


P X2 


p-^ p^ 

where the second terms on the right hand side constitute the general solution of the 
homogeneous problem. Now we substitute this into the limit of the conservation 
equation 


to obtain 


dt{pf - Pg) + v{f - ^)dn^ = 0 , 

dt (p- 


= V. 


Dip + ) 'Vxp 


where we have introduced the positive dehnite symmetric matrices 


Di = V vxi dn and D2 = / n <8 VX2 dn . 


This can be written as the nonlinear diffusion equation 

dtm = Vx ■ {D{m)Vxm) , 
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for the new unknown 

m = p -, i.e. p{m) = - + \/ m? + 4 ^ , 

where we have introduced the diffusion matrix 

Dim) = (Oipimf + x- 7 -^-• 

\ pim)-^J 2pim)—m 

The unknown m is the zeroth order approximation of the difference of the position 
densities of / and g. 


2. Long time properties 

In this section we study decay to equilibrium for solutions of o and of the 
linearized problem ([7]), ([ 8 ]). In order to estimate the decay towards the equilibrium 
quantitatively we introduce the scalar product, weighted with the steady state 
measure, 

+ili^']dvdx, withF,= f* 


( 10 ) 


(Fi,F 2 ) = 



't3 Jr^ \PooXi X2 

Throughout this article we denote by || • || the norm induced by this scalar product. 
The orthogonal projection onto the null space of the linearized collision operator L, 


( 11 ) 


IF = 


PooXlPg 

~X2pf-poog^ 
that is the space of local equilibria, is given by 


with F = 


( 12 ) 


HF := 


Pf- Pg / PooXi 


pL + 1 


-X 2 


Straightforward computations yield 
-{IF,F) = 


IT 3 


(/ - PfXl? ^{ 9 - PgX2)^P\ 


2 .,2 

00 


XiPIo 


X2 


dn 


Pf 

— 1“ PgPa 

Poo 


dx , 


and 


||(i-n)F|r 



(/ - PfXlY ^ {9 - PgX2)^Poo \ ^ Poo 


XlPc 


X 2 


pL + 1 VP< 


Pf 


T PgPa 


dx , 


implying the microscopic coercivity estimate 


(PI, H3) 


- {IF,F) > min{poo, 1/Poo} ||(1 - n)F|l^ 
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This gives a quantitative estimate of the decay towards the local equilibrium intro¬ 
duced by the linearized collision operator. In the spatially homogeneous situation 
such an estimate is enough to prove exponential decay to equilibrium for the lin¬ 
earized equation. 

For spatially non homogeneous situations we expect the densities to become con¬ 
stant as these are the only local equilibria that also annihilate the transport part of 
the equation. The complete relaxation mechanism can be seen as a combination of 
local relaxation in the velocity direction by the collision operator and an interplay 
between mixing by the transport operator and confinement in our bounded spatial 
domain. In the following we will study two recent methods mm and 113 ) to esti¬ 
mate the decay in the spatially non-homogeneous situation. Both rely on properties 
of the linearized collision operator that we will verify in the sequel. Concerning 
their numbering, (Hn) and (Pm) in this work correspond to (Hn) in [12] and, re¬ 
spectively, (Hm) in [8]. The microscopic coercivity property ( |P1, H3 ) is needed in 
both approaches. 


2.1. Coercivity in weighted H^. 

When studying coercivity of the collision operator we saw that the operator provides 
coercivity only with respect to the velocity distribution. The strategy in [12] is to 
transfer some of this dissipation effect in the velocity to the spacial variable by using 
the mixing effect of the transport operator. This method was mainly inspired by 
discussions with C. Villani (see [15]) and results by Y. Guo (see for example [10] or 
the references in Ha)- 

We start by writing the linearized equation ([7]) in the abstract form 
(13) f + = 

with F = (/, g)'^, with the transport operator 


TF = 



(v • V./\ 
\v ■ VccgJ ’ 


and with the linearized collision operator L given in flTTl) . The regularizing effect of 
the transport operator can be quantified by looking at the time evolution of mixed 
derivative terms 


(14) ;|(V,F, V,F) = -(V,F, V,TF) + 2(V,LF, V,F) = -\\V^Ff + rest, 

where the derivatives are to be understood in the sense VF = (V/, Vg)'^. Thus the 
transport operator provides some damping in x-derivatives provided the rest terms 
can be compensated by boundedness and coercivity in the velocity direction of the 
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collision operator. We introduce the Sobolev-spaces with weights Xi and X 2 
and the norm 

11^11^1 = ll^ll' + l|V.F||^ + ||V,F|^ 

In pS] a norm ||F||^i = ||F||^i + 6{VxF, VyF) is constructed that is, for 6 small 
enough, equivalent to the hT^-norm but uses these mixed derivatives to get a co- 
ercivity estimate also for the spatially non-homogeneous situation. Exponential 
convergence in is then derived from 

( 15 ) illFlg. <-rllFfn,, 

and convergence in with the correct weight follows from the equivalence of the 
norms. Since the essential mixing effect of the transport operator is quantihed by 
means of fll4p this approach is suitable only for proving coercivity in a space of higher 
differentiability, more precisely spaces, weighted with the the steady state (or 
Gibbs-) measure. This is a restriction of the method on one hand but on the other 
hand the method can also be employed to prove convergence in higher derivatives. 
This allows to use embedding theorems and lends itself to studying nonlinear equa¬ 
tions in the perturbative regime. In the model we study here, however, we have 
good a priory bounds already and thus there is no need to go beyond as we will 
see in the proof of Theorem 01 


To control remainder terms in the time evolution of the "H^-norm structural as¬ 
sumptions on the linearized collision operator are made. We will verify these in 
the sequel. For a collision operator that is local in x, as is the case for our model, 
the spatial derivatives cannot increase as can be seen by using ( PI, H3 ) on the 
derivatives. To get decay of the velocity derivatives it is essential that the reaction 
operators can be split L = K — A, into a “loss” part 

7 ^ 

\9 / \Poo9 / 

that is (in our case trivially) coercive on the L^-level 


A 


—/ 

Poo ■' 


(Hl/2) (Af’,F)>min{p„„,jL}||F||y 

and a gain part 

\3J \-^X2PfJ 

which is regularizing in v as long as xi X 2 are regular. Indeed, for V„ acting 
component-wise in the two functions, a standard interpolation argument yields 

(H2) V<5>0:3G: \{V,KF,V,F)\ < S \\V,Ff + C \\Ff , 
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where C depends on ||V^XilU°° and || 

Since the method relies on type estimates actnally coercivity in of the loss 


part is necessary. In onr case the same estimate as in (Hl/2) resnlts in 


(Hl/3) 


{V„AF. V„F) > min|p„, l|V„F|r 


where in general negative terms of lower order derivatives are allowed bnt not needed 
in this case. 

Control of terms of the type (VLF, VF) is ensnred by the fact that the linearized 
collision operator is bonnded in the sense that 


(Hl/1) 




PgU + 


fu 

PIoXI 


+ PfV + 


PVplo 

X2 


dn dx 


< 4max{poo, 




This can be seen by nsing J Xidn = J X2dv = 1 and Canchy Schwarz in estimates 
of the type 



Properties (Hl/1) to (Hl/3) together with (IH2D . ( |P1, H3D ensnre that we can 
nse the main theorem from [ 12 ] to derive the following convergence resnlt for the 
linearized problem; 


Theorem 2. Let xi o,nd X 2 be in (M^) and the initial data /o be in H^{dv/xi) 
and go in H^{dv/x 2 )- Then the solutions f,g of the linearized problem ([^ subject 
to initial conditions f{t = 0) = fo, g{t = 0) = go exist globally and converges 
exponentially to the equilibrium distribution. For fjsfjg^sifo — go)dvdx = 0 the 
equilibrium is zero and we have 

ll/(•,•,^)II^Il(d,;/xl) + ll^(•,•,^)II^Il(d^/x2) < C'exp(-rf), 

where the rate r depends on the constants in the estimates tagged with (H1)-(H3) 
and the constant C also depends on the norm of the initial data in and of Xi 

y^l.oo 































10 


LUKAS NEUMANN, CHRISTIAN SCHMEISER 


Remark 3. Convergence in higher order Sobolev-spaces can he derived straightfor¬ 
wardly provided the estimates (Hl/3) and flH2p can be generalized to higher order 
derivatives, as is easily verified to be the case for our model. 


This feature is usefull mainly in applying the results to the nonlinear system in a 
perturbative setting. Control of the bilinear contribution in the interaction is given 
by applying chain rule, Holder inequality and using Sobolev-embedding to lower the 
exponents in the norm to two again. In Dimension 3 we see that 

(H4) yk>2: <C 

\Pf9j \9J 


Now the exponential decay of the H^- norm in the linearized situation can be used 
to conclude that for initial data close to the global equilibrium in we have 
convergence to the stationary state. 

Here however we want to give a stronger result - in the sense that less regularity 
is necessary - by using the a priori bounds of Theorem [T] 


Theorem 4. Let Xi o-nd X 2 be in (M^) and the initial data /o be in H^{dv/xi) 
and go in H^{dv/x 2 )- Moreover let fo and go satisfy the assumptions of TheoremU\ 
with 7 i and 72 small enough. Then the solution to equations o with initial data fo, 
go exists globally in time and converges to the unique stationary state exponentially, 
more precisely 

ll/(•) •) ~ PooXi\\h^(Av/xi) + II5'(l •) 1) — X 2 / Poo\\h^(Av/x2) <Ce 

where the constants r depends only on the estimates (H1)-(H3) and C depends also 
on 7 i, 72 as well as the -norms of Xi- 


Proof. We use the results for the linearized collision operator and regard the differ¬ 
ence to the nonlinear one as a small perturbation. This difference is given by 

'ipg - l/Poo)(/ -PooXi)' 
ipf Poo){,9 X2/ Poo) 

and, by Theorem (H pf is close to poo and Pg close to 1/poo, yielding 

(16) ||Q(/,p)-LF||<7||F|| 


Q(/,p)-LF = - 


with a small constant 7 depending on 71 and 72. We set Fao = (PooXi; X2/P00) and 
apply our results for the linearized collision operator to F — Foo, where F solves ([1]). 
Note that since L{Foo) = 0 we have for evolution under the full equation ([1]) 


dt (F - Foo) + T (F - Foo) = L (F - Foo) + Q{f, p) - LF . 

Thus equation flTS]) becomes, including the derivative of the nonlinear interaction, 
(17) ;|||F - Fooll^i < -r||F - Fooll^i + 2||F - FooH^i||Q(/,p) - L{F)\\ni . 
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Estimate ffT^ holds for velocity-derivatives straightforwardly. For spatial derivatives 
we use the bounds of Theorem [T] and the multiplicative structure of the nonlinearity 
to derive 


V, {Q{f,g)-L{F))\\ < 

(/-PooXl)\^ f{pg 

X 2 /P 00 )) ^ 


{9 


1/ Pot 

(Pf ~ doo) 


+ 


(Pg 1/ Pot 

iPf Poo) 


if - PooXl) 
{9 - X 2 /P 0 O) 


< 


7 


^ (Xl{pg-^lpoo)\ 

1 

^ f if - PooXl) \ 

""V X2{Pf-Poo) J 

+ 7 

\i9-X2/Poo)) 


< 


\F-Fr 


coW'H} 1 


where 7 is a small constant that changes form line to line and we used the fact that 
J Xidn = 1 = J X 2 dn togehter with Chauchy-Schwarz in the last estimate. We infer 

\\Q{f,g)-L{F)Ur<^\\F-F^Ur 


and using this in flTTD yields the exponential convergence as long as 71 and 72 , and 
thus 7 , are small enough. □ 


2.2. Coercivity in a weighted L^-space. In this section we apply the abstract 
convergence theory of [ 8 ]. This approach does not use derivatives to quantify the 
mixing effect of the transport but rather uses a modihed entropy functional. We 
again start with the linearized equation ([7]) in the abstract form 


dF^ 


LF. 


It is easily seen that L is symmetric and T is skew-symmetric in the Hilbert space 
■H, defined as the weighted L^-space with the scalar product given in flTOl) . In this 
space the map H, defined in ffT 2 l) . is the orthogonal projection to the null space of 
L. 

The approach of [ 8 ] relies on the modified entropy functional 

H[F] := + (5(AF, F ), with A := (1 + (Tn)*Tn)-FTn)* , 

and with a small positive constant 5. The time derivative of the modified entropy 
along solutions of flT5]) is 

(18) = (LF,F)-(5(ATnF,F)-(j(AT(l-n)F,F) + (5(TAF,F)+(5(ALF,F). 

CJ .6 


The hrst term on the right hand side suggests that the microscopic coercivity esti¬ 
mate (PI, H3) is one of the necessary ingredients. Since the operator ATH can be 
interpreted as the application of the map z i-7 to (Tn)*Tn, the second condition 
called macroscopic coercivity is also plausible: There exists \m > 0, such that 


(P2) 


l|TnF||">AMl|nFf. 
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As a consequence of (PI) and (P2), the sum of the hrst two terms in the entropy 
dissipation above is coercive. It controls the remaining three terms, if the operators 
appearing there are bounded and act only on the microscopic part (1 — n)F of the 
distribution. For the operator A this is guaranteed by the algebraic condition 

(P3) nrn = o, 

called parabolic macroscopic dynamics (see Section 1.3). The hnal condition is the 
boundedness of the auxiliary operators: There exists Cm > 0, such that 

(P4) ||AT(1 - n)F|| + IIALFII < Cm\\{1 - n)F|| . 

Boundedness results of the same form for A and TA hold as a consequence of (P3) 
(see Lemma 1 of [H]). The former leads to coercivity of H[F] for 6 small enough. 
We formulate the abstract convergence result from [8]: 

Theorem 5. Let fPli- flP4p hold and let F be a solution of flT^ subject to F{0) = 
Fq G TL. Then there exist constants X,C > 0, such that 

\\F{t)\\<Ce-^^\\F4. 

For our problem it remains to verify (P2)-(P4). A straightforward calculation 
shows that (P2) is equivalent to 



WxU^DqV xUpdx > Am 



M^dx, 


with Up = Pf ~Pg satisfying updx = 0 by ([8]). Thus, by the positive dehniteness 
of Do = (p^ + + D 2 ), (P2) is a consequence of the Poincare inequality 

on T^. Since 


V ■ VxUp f Pcx)XA 

pL +1 V -X2; ’ 


and since the application of 11 involves an integration with respect to n, the as¬ 
sumptions ([3]) imply flP3p . By (Hl/1) the linearized collision operator L is bounded. 
For the verihcation of (P4) it is thus sufficient to prove the boundedness of AT or, 
equivalently, of its adjoint (AT)* = —TA* = —T^n(l -|- (Tn)*Tn)~^. The equation 
G = (1 + (Tn)*Tn)-iF implies 


Ug — Vx ■ {DqVxUg) — Up . 

The norm ||T^nG|| is equivalent to the L^(T^)-norm of whose boundedness 

in terms of the L^(T^)-norm of up (and therefore in terms of ||T||) is a consequence 
of elliptic regularity. This proves (P4) and completes the proof of exponential decay 
to equilibrium for the linearized problem. 
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Theorem 6 . Let Q hold, let /o G L‘^{dv/xi), go G L'^{dv/x 2 ), CLnd let J ^3 fg^sifo — 
go)dvdx = 0. Then the solution of ([7]) subject to fit = 0) = fo, g(t = 0) = go, 
satisfies 

ll/(■) ■) + ||5'(') ■) ||L2(dli/x2) — 5 

with positive constants C and A. 

Since the maximum principle estimates of Theorem [1] already imply a stability 
(but not asymptotic stability) result for the nonlinear problem, the decay result can 
be extended to a local result for the nonlinear case by the same method. 

Theorem 7. Let (jS]) hold and let fo and go satisfy the assumptions of Theorem [I] 
with 7 i and 72 small enough. Then the solution of the initial value problem O, (0 
satisfies 

PooXl II L'^{dv/xi) ~h ||l?(') ■)^) X 2 / Poo 11’ 

with positive constants C and A. 

Proof. We start by writing the problem in terms of the unknown F = (/ — PooXi) S' ~ 
X 2 /Poo)^- Then we proceed as above producing the entropy decay relation (ITH]) . 
however with LF replaced by 

W.9)= 

VX 2 - pfgj 

Since, by Theorem [H p/ is close to poo and Pg to 1/poo we have again (cf. 

||Q(/,p)-LF|| <7||F||, 

with a small constant 7. Therefore the entropy dissipation for the nonlinear equation 
is a small perturbation of the entropy dissipation of the linearized problem, which 
does not destroy its coercivity. □ 

3. Rigorous macroscopic limit 

Our next goal is to validate the macroscopic limit carried out formally in Section 
0 We start from the rescaled system ([2]), subject to initial conditions, where the 
data satisfy the assumptions of Theorem [T] The results of Theorem [T] remain valid 
with the e-independent L°°-bounds. We start by exploiting entropy decay. 

Lemma 8 . Let the assumptions of TheoremUlhold, let {f,g) be the solution of (P|), 
on]) fore > 0 , and define the micro-macro decompositions f{x,v,t) = pf{x,t)xi{v)-\- 
ef-^{x,v,t), g{x,v,t) = pg{x,t)x 2 {v) + ^g^{x,v,t). Then f^ and g^ are bounded 
uniformly in e in x x ( 0 ,00), dx dn df/vO, ? = 1,2, respectively, and 

= 0{e) m L2(T3 x ( 0 , 00), dxdf). 
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Proof. The proof is based on the entropy decay relation 
e^dH{f,g) 


2 dt 



/T3 JR3 


X1X2 fl - In ^^dn'dnda;. 

V X 1 X 2 / X 1 X 2 


Since the entropy is nniformly bonnded in e and f, nsing {^/a — 1)^ < — 1) In a, 

we derive 



fg' 


10 JT 3 JR 3 


X 1 X 2 \ - T - 1 dv'dvdxdt = O(c^) 


X 1 X 2 


Using the micro-macro decomposition and expanding the sqnare we hnd 



r- -TT--— ^ ea eb e^ia — hY 

\/ (1 T ecz)(l -|- e6) — 1 -|- — -|- — —-^- 

2 2 4(y/(Y+7a)(r+7^-h 1ea/2e6/2) 

with a = f^/{pfXi), b = g^'/{pgX 2 )- Since 1 + ea = //(p/Xi), I + eb = g'/ipgXY^ 
the estimates from Theorem [T] can be nsed to obtain 1 -|- ea, 1 -|- e6 < £^£±22 -with 
the conseqnence 

Using this in flT^ . we obtain 

e^(Poo - 7i) 


0{e) = I{t)> {,/pJpf-l)dx+ . 

Jt3 4(poo + 72 )'^ Jt^ 

completing the proof. 



^ ^ dn da;, 

Xl X 2 / 


□ 


With this basis we now follow the procednre of [H]. In particnlar, we nse an 
averaging lemma, which can be proved similarly to Lemma 3.2 in [13]. We give the 
short proof for completeness. 

Lemma 9. Let x satisfy ([3]), let f and h lie in subsets o/L^(T^xM^xM, da; dn dt/x{v)), 
uniformly bounded in terms of the small parameter e, and let edtf + v ■ Vxf = h. 
Then pf is bounded uniformly in e in L^(M; 
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Proof. We represent / by the Fourier transform with respect to t and by the Fourier 
series with respect to x: 

f{x,v,t)= 

^eT3* 


with the lattice T^* dual to the torus T^, implying 

zf = —ih , with z = er + V ■ f . 

For each A > 0, we introduce a smooth, nonnegative real function fjxlz) < 1, 
satisfying 'ipxi^z) = 0 for \z\ < A and 'ipxi^) = 1 for 1^1 ^ 2A. Now we estimate, using 


|P/I < 


Hav 


+ 


(1 -'ilJx)fdv 


< 


ll^l|L2(d^/;^) + ( / {l-fjxfxdv) ||/||L2(d„/^) 

R3 2 ; y \jj^3 y 

— ^ll^l|L2(di;/x) fo \]^/ ll^llL2(di;/x) • 

With the optimal choice A = obtain 

completing the proof. 


\ 1/2 


□ 


Theorem 10. Let the assumptions of TheoremU\hold. Then as e —)■ 0 the solution 
if, 9 ) 0 / (I2D, ® converges to {pXi,X 2 /p) in Ll^{T^ x x (0, cx)), dx dn df/xi) x 
X X (0, cxo), dx dn df/x 2 ); when restricting to subsequences, where p G 
L°°(T^ X (0, cxo)) satisfies poo — 7i < P < Poo + 72- Furthermore there exist Ji, J 2 G 
L^(T^ X (0, oo) fi such that 

dt (p -) + Vx • (— </2) = 0 , -Ji = —iAiVxP, pJ2 = —D2^x—, 

\ pj P P 

hold in the sense of distributions. 


Proof. Because of the boundedness of pf and Pg and of Lemma [H / and the function 

h := -= Xi(l + x/pjp^) -^- Pgf 

satisfy the assumptions of Lemma [9] with x = Xi (after even extension to t < 
0). With an analogous argument for g we obtain Pf,Pg G L^((0, cxo); 
uniformly in e. 










16 


LUKAS NEUMANN, CHRISTIAN SCHMEISER 


The conservation law 

( 20 ) dt{pf - pg) + v{f^ - = 0 , 


the observation 


( 21 ) 



9^)dv 


< “s/tiDi 11/ + v^trZ^ ll^f ||-L2(dii/x2)) 


and Lemma [H] imply pf — pg ^ oo)] which, after interpolation with 

the averaging resnlt gives 

e „ 

Pf — Pg E iL2(i+s) ((0, oo) X M ) nniformly in e . 

As a conseqnence, for each 0 < a < 6 and compact K C a snbseqnence oi pf — pg 
converges strongly in L^((a, 6 ) x iL) as e —)■ 0. Since the same is trne for ^Pfpg —)■ 1, 
it also holds for pf and pg individnally as a conseqnence of the L°° bonnds. Another 
application of Lemma [ 8 ] completes the proof of the convergence statement. 

For the derivation of the limiting problem, we pass to the limit in (120]) in the 
distribntional sense, denoting the weak limits of vf^dv and vg^dv, which 
exist becanse of fl2T]) . by A and J 2 , respectively. Now we mnltiply the eqnation for 
/ by n/e and integrate with respect to v obtaining 

edt / vf^dv + V 3 ; • / V ® vfdv = —pg / vf^dv . 

JM? Jr^ 

By the nniform-in-e bonndedness of Jj^g v <8 vfdv (conseqnence of Theorem [1] and 
([3])) and by the strong convergence of pg, we can pass to the limit, leading to the 
desired eqnation for Ji. For J 2 we proceed analogonsly. □ 


References 

[1] Bisi, M., Desvillettes, L., From reactive Boltzmann equations to reaction-diffusion sys¬ 
tems, J. Stat. Phys. 124, no. 2-4 (2006), pp. 88I-9I2. 

[2] Bothe, D., Hilhorst, D., a reaction-diffusion system with fast reversible reaction, J. Math. 
Anal. Appl. 286 (2003), pp. 125-135. 

[3] Carrillo, J., Desvillettes, L., Fellner, K., Fast-reaction limit for the inhomogeneous 
Aizenman-Bak model. Kinetic and Related Models 1, no.l (2008), pp. 127-137. 

[4] Carrillo, J., Desvillettes, L., Fellner, K., Rigorous derivation of a nonlinear diffusion 
equation as fast-reaction limit of a continuous coagulation-fragmentation model with diffusion. 
Comm. Part. Diff. Eq. 34 (2009), pp. 1338-1351. 

[5] Choquet, L, Degond, P., Schmeiser, C., Energy-transport models for charge carriers 
involving impact ionization in semiconductors. Transport Theory and Statistical Physics 32 
(2003), pp. 99-132. 








KINETIC REACTION MODEL 


17 


[6] Degond, P., Nouri, a., Schmeiser, C., Macroscopic models for the ionization in the 
presence of strong electric fields, Transport Theory and Stat. Phys. 29 (2000), pp. 55I-56I. 

[7] Dolbeault, J., Mouhot, C., Schmeiser, C., Hypocoercivity for kinetic equations with 
linear relaxation terms, C.R. Acad. Sci. Paris 347 (2009), pp. 5II-5I6. 

[8] Dolbeault J., Mouhot C., Schmeiser C., Hypocoercivity for linear kinetic equa¬ 
tions conserving mass, Trans. AMS, electronically published on February 3, 2015, DOI: 
http://dx.doi.Org/10.I090/S0002-9947-20I5-060I2-7 (to appear in print). 

[9] Golse, F., From kinetic to macroscopic models, in Kinetic Equations and Asymptotic Theory, 
B. Perthame, L. Desvillettes (Eds.), Series in Appl. Math. 4 (Gauthier, Villars), 2000, pp. 41-- 
126. 

[10] Guo, Y. The Vlasov-Poisson-Boltzmann System near Maxwellians. Comm. Pure Appl. Math. 
55 (2002), 1104-1135. 

[11] Hilhorst, D., van der Hout, R., Peletier, L.A., Nonlinear diffusion in the presence of 
fast reaction. Nonlinear Anal.: Theory, Meth. & Appl. 41, no. 5-6 (2000), pp. 803-823. 

[12] Mouhot, C., Neumann, L., Quantitative perturbative study of convergence to equilibrium 
for collisional kinetic models in the torus. Nonlinearity 19, No. 4 (2006), pp. 969-998. 

[13] POLEWCZAK, J., The kinetic theory of simple reacting spheres: I. Global existence result in 
a dilute-gas case, J. Stat. Phys. 100 (2000), pp. 327-362. 

[14] POUPAUD, F., Schmeiser, G., Charge transport in semiconductors with degeneracy effects. 
Math. Meth. in the Appl. Sci. 14 (1991), pp. 301-318. 

[15] ViLLANi, C., Hypocoercivity, Memoirs of the AMS 950, 2009. 


L. Neumann 

Universitat Innsbruck 
Technikerstrasse 13 
6020 Innsbruck 
Austria 

E-MAIL: Lukas.Neumann@uibk.ac.at 


C. Schmeiser 

Universitat Wien 
Oskar-Morgenstern-Platz 1 
1090 Wien 
Austria 

E-MAIL: Christian.Schmeiser@univie.ac.at 



